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Abstract - We study the repulsive Hubbard model on an anisotropic honeycomb lattice within 
a mean-field and a slave-rotor treatment. In addition to the known semi-metallic and band- 
insulating phases, obtained for very weak interactions, and the anti-ferromagnetic phase at large 
couplings, various insulating spin-liquid phases develop at intermediate couplings. Whereas some 
of these spin liquids have gapless spinon excitations, a gapped one occupies a large region of the 
phase diagram and becomes the predominant phase for large hopping anisotropics. This phase 
can be understood in terms of weakly-coupled strongly dimerized states. 
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One of the most salient features of graphene, an atom- 



^ically thin graphite sheet with carbon atoms arranged in 
'— !a honeycomb lattice, is certainly the relativistic charac- 
^.^ter of its low-energy electronic excitations [l]. An im- 
^portant issue of present research efforts is a deeper un- 
derstanding of the role of interactions in this novel two- 
dimensional material and of the persistance of the semi- 
'\! metallic (SM) phase when such interactions become rele- 
^'.vant. As a consequence of the poor screening properties 
^>>Oof electrons in the SM phase, the Coulomb interaction re- 
mains long-range and may be described in terms of an ef- 
fective fine-structure constant that turns out to be density- 
• •independent, ao = e^ /hev-p ~ 2.2/e (for recent reviews on 
.^electronic interactions in graphene, see Refs. [2] and [sj). 
^ Here v-p ~ c/300 is the Fermi velocity in terms of the speed 
$-H of light c, and e is the dielectric function of the surround- 
^ ing medium. Whereas lattice-gauge theories predict a flow 
to strong coupling above some critical value ~ 1.1 4|, 
even suspended graphene seems to be weakly correlated, 
with a stable SM phase [S] . A reason for this effective flow 
to weak coupling may be an intrinsic dielectric constant in 
graphene due to virtual interband excitations f6]. Indeed, 
recent renormalization-group studies confirm this picture 
of weakly-interacting electrons in graphene [t] . 

A perhaps more promising system for the study of the 



interplay between strong (short-range) correlations and 
the relativistic character of Dirac fermions may well be 
a gas of cold fermionic atoms trapped in an optical hon- 
eycomb potential [8j[9]. As compared to graphene, such a 
system has several advantages. First, neutral atoms resid- 
ing on the same site exhibit short-range interactions which 
can often be tuned over orders of magnitude and turned 
repulsive or attractive by using Feshbach resonances 10 



Indeed, recent experiments have proven the feasibility of 
implementing optical honeycomb lattices and probing in- 
teraction physics in the context of bosonic ^^Rb atoms 
[TT] . Second, the hopping rates between neighboring sites 
are easily controlled by the laser configuration and beam 
intensities as both determine the depth and position of the 
optical potential wells. Rather moderate changes in the 
laser intensities can significantly imbalance the tunneling 
rates, leading to situations ranging from weakly-coupled 
zig-zag linear chains to weakly-coupled dimers [9]. This 
situation needs to be contrasted to graphene, where un- 
physically large lattice distortions are required to obtain 
the limit t' ~ 2t [12] , where novel physical phenomena are 
expected [l3] . 

Starting from the repulsive fermionic Hubbard model 
(RFHM) for spin- 1/2 particles with onsite interactions 
[/>0 and identical nearest-neighbor hopping rates t, mean- 
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field calculations for a half-filled lattice and zero temper- 
ature predict 

a SM phase at small values of U jt and an anti- 
ferromagnetic (AF) phase developing above — 2.2t 
[l4^jl5j. Quantum Monte Carlo (QMC) calculations con- 
firm this SM-AF transition but at a higher transition point 
Uc ~ 4.5t 14 16 while dynamical-mean-field estimates 
yield even larger values, Uc ~ lOt 17 . Most saliently, 
recent QMC investigations have revealed an intriguing in- 
sulating spin-liquid (SL) phase, with localized charges but 
no spin ordering, that emerges below the AF transition 
and that may be related to the exotic Mott insulator 

The transition 



16 



identified in slave-rotor studies [18f 20 
points derived from the slave-rotor theory also occur 
at globally smaller values of U jt than in QMC calcula- 
tions. 

Here, we study the RFHM on a half-filled honeycomb 
lattice where one of the nearest-neighbor hopping param- 
eters t' is larger than the other two t. Such a lattice re- 
veals an astonishingly rich phase diagram (see Fig. [T]), in 
which SL phases predominate at large values of t' jt. At 
?7 = 0, this system develops a topological phase transition 
between a semi-metal and a band insulator (BI). Indeed, 
at t' — 2t, the two Dirac points responsible for the SM 
phase merge and eventually disappear while a band gap 
These topological properties are prominently 



opens 

refiected in the SL phase, where the spin excitations ac- 
quire a gap as a function of renormalized hopping parame- 
ters. Finally, we show that the SL phase finds a compelling 
interpretation in terms of weakly-coupled dimer states in 
the limit > i. 

As the honeycomb lattice is made of two shifted identi- 
cal triangular sublattices A and B, the kinetic term of the 
RFHM, H = Ho + V, reads 



Hq 



•J 



{i',j')cr 



where a'^'J and b^^J annihilate (create) a fermion with 
spin a =t; i on the sites of the A and B sublattices, re- 
spectively. It also distinguishes the hopping amplitude 
t' > 0, chosen along neighboring sites linked by the vector 
Ci — ae^, which differs from the other two t > chosen 
along neighboring sites hnked by C2 = a(— e2./2 + \/3e.y/2) 



or by C3 = -a(ej;/2 
stant [see Fig. [2]ja)]. 
term reads 

V = 



y/3ey/2), a being the lattice con- 
At half-filling, the onsite interaction 

U 



E(^ 



■it ■ 



HI 



(2) 



with U > 0, the summation running over both sublattices 
and fiia = Oi^cLicr or bl^bia being the corresponding num- 
ber operators. Here we consider a balanced population 
between spin-f and spin-|. When [/ = 0, the Hamilto- 
nian is readily diagonalized in reciprocal space, and one 
obtains the two energy bands = ilT^^lj in terms of 
the weighted sum of phase factors 7*° = tg'' + t'g''^ with 
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Fig. 1: (Color online) Mean-field phase diagram obtained for 
the anisotropic repulsive Hubbard model on the honeycomb 
lattice at zero temperature and half-filling. The upper thin- 
ner blue line delineates the AF phase, whereas the two linear 
functions (the dashed lines) for [7c are obtained from simple 
scaling arguments (see text). The lower thick green line indi- 
cates a second-order phase transition between the SM and the 
gapless SL of the third solution (GLSL), whereas the dotted 
line are first-order transitions, which consists of three sections: 
the black round dots for LOOS < i! jt < 1.02 between the GLSL 
of the third solution and that of the chain solution, the black 
squares for 1.02 < t' /t < 1.31 between the SM and the GLSL 
of the chain solution and, above t' /t ~ 1.31, the red triangles 
between the SM (or BI) and a gapped SL (GDSL) that consists 
of decoupled dimers on the bonds with t' . The curve across the 
black squares and red triangles is drawn to direct the eyes. 



recovers the usual Dirac points located at K and K', as 
depicted in Fig. ^h). The Dirac points move towards 
the point M as t' is increased, where they finally merge 
2t (8)[22 



+ e 



and 



_ ^ik ci 



When t' ^ t, one 



At half-filling, this topological phase 
transition separates a SM phase (for t' < 2t), with two 
Dirac points, from a BI (for t' > 2t), the insulating gap 
Al = t' — 2t opening at the point M. 

For this half-filled lattice, one expects a Mott insulat- 
ing (MI) phase at large U in which each lattice site is 
exactly occupied by one atom. As a consequence of the 
superexchange interaction At'^/U, the spins of the parti- 
cles are further frozen in an AF Neel state that breaks the 
lattice inversion symmetry, with e.g. spin-f atoms on the 
A sublattice and spin-j, on the B sublattice. This state 
is described by the staggered mean-field order parameter, 
{hi^ — fii^) = A if I G A, and — A if i G -B. The mean- 
field approximation fii^nii w [(1 ± A)ni-f + (1 ^ A)ni^]/2 
in the interaction term V leads to a quadratic Hamilto- 
nian which is diagonalized through a Bogoliubov-Valatin 
transformation. The corresponding excitation spectrum 
is = zb^A^j -t- |7fc|2. The Mott gap is found to be 
Am — UA /2 and is therefore intimitely related to the AF 
order parameter, i.e. it opens as soon as the AF order sets 
in. In the thermodynamic limit, the self-consistent gap 
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Fig. 2: (a): The honeycomb lattice and its diamond-shaped 
unit cell delineated by the Bravais lattice vectors ai^2- Black 
(resp. grey) sites belong to sublattice A (resp. B), a being the 
distance between nearest-neighboring sites. The hopping am- 
phtude t' connects sites linked by the vector ci while t connects 
sites linked by 02,3. (b): The reciprocal lattice is spanned by 
the primitive vectors bi,2, and has the diamond-shaped prim- 
itive unit cell as shown. The F point is the center of the unit 
cell. In the isotropic hmit, the two inequivalent Dirac points 
are located at K and K'. 



equation at zero temperature then reads 
d^k U 



FBZ ^ v/C^2A2+4|^fe|2 



(3) 



where n = Stt'^ / {Sy/Sa'^) is the area of the first Brillouin 
zone (FBZ). The transition line Uc delineating the AF 
phase (blue solid line in Fig. [l]) is obtained from Eq. ([3| 
by setting A = 0. Notice that the value Uc — 2.2t, ob- 
tained at t' = t, agrees with former mean-field calculations 
14 15 . The critical value Uc is however shifted to larger 
values when t' increases. This shift may be understood 
from a simple scaling argument when considering the full 
band width W — 2{2t + t') of the non-interacting case. 
Since all particles are localized in the AF phase, one needs 
to consider states with energies up to ~ W . The critical 
value for the isotropic case t' = t may now be expressed as 
Uc/W ~ 0.37. If we consider this value to remain constant 
when varying t' , one obtains Uc — 0.74t' -I- 1.48i, a rela- 
tion that describes to great accuracy the transition line 
for t' < 2t (Fig. [1]). At larger values of t' , the mean-field 
Mott gap Am competes with the band insulator gap Ai, 
and the transition line may be understood qualitatively as 
the line where both gaps are equal. This yields the asymp- 
totic behavior Uc = 2{t' — 2t) ^ 2t' in the large t'-limit 
(dashed line in Fig. [I]). 

In order to decouple the Mott transition from the AF 
phase and to study possible intermediate phases, more so- 
phisticated methods than a simple mean-field treatment 
are required. Apart from QMC calculations 16 , an in- 
termediate MI spin-liquid phase has been identified in the 



honeycomb lattice within a slave-rotor treatment 18 ■ 20 
where the fermion operators aj^ and bja- are viewed as 



products of two auxiliary degrees of freedom, exp{i9j)fjcr- 
We adopt, here, a U(l) slave-rotor treatment that is ex- 
pected to provide qualitatively correct results within the 
mean-field approximation. However, if one aims at an 



effective low-energy theory for the spinon degrees of free- 
dom, one needs to take into account the coupling to an 
SU(2) gauge field as discussed in Ref. 19 . The bosonic 
"rotor" field 6j is conjugate to the total charge at site j, 
described by the angular momentum Lj = idg^ , and the 
fermion operator fj^ carries the spin ("spinon"). As this 
procedure artificially enlarges the Hilbert space, double 
counting needs to be cured by imposing the constraint 



(4) 



at each site j. Whereas the rotor and spinon fields are 
coupled via the kinetic Hamiltonian T-Lq, the interaction 
term is described solely in terms of the rotor degrees of 
freedom, V = iU/2)J2j L'j, the particle-hole-symmetric 
form of which also renders the chemical potential fi — 0. 
Following Ref. [20j, the rotor and spinon degrees of free- 
dom may now be decoupled within a mean-field treat 
ment by defining the averages Qg'' 
Q^p = (exp(-i%)), where % = 9i 
nearest neighbors connected by Ci for the primed aver- 
ages and by 02,3 otherwise. Notice that the mean-field 
parameters along C2.3 are assumed to be equal, hence the 
corresponding symmetries are not broken here. The de- 
coupled mean-field Hamiltonian may then be written as 
H ~ He+Hf, with 



HAflha) and 
'j and i and j are 



t' g^e-'«-'.' + (h.c.) 
2 



(5a) 



2 ^ V ' U 

i ' 

-Y.^Jlh., (5b) 



where hi is a local Lagrange multiplier ensuring the con- 
straint Q . At half-filling, particle-hole symmetry imposes 
that hi = Q. The Mott transition may now be inter- 
preted in terms of rotor condensation. Indeed, in the 
rotor-condensed phase, the phase 6j is fixed and the num- 
ber of particles (or the angular momentum Lj) therefore 
fluctuates on the lattice sites. This corresponds to the 
SM phase for t' < 2t or the BI for t' > 2t. In the MI 
phase, however, there is no spin ordering since the spinon 
Hamiltonian % / has no interaction term. 

Most salicntly, "Ho can be mapped onto T-Lf through 
[f ,t) — ^ {t' Q'j: , tQ f) , such that the low-energy spinon ex- 
citations are governed by the same topological properties 
as the non-interacting Dirac fermions. In particular, sim- 
ilarly to 7*", we find the weighted sums of phase factors 



and 



(6) 



(7) 
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Fig. 3: (Color online) Spinon dispersion obtained from mean- 
field Hamiltonian T-Lf. (a) For t'Q'^ = tQf, the dispersion 
consists of two distinct Dirac points, and one obtains a gapless 
spin liquid, (b) For t'Q'^ — 2.5tQf, the dispersion relation is 
gapped for any vector k (gapped spin liquid). 



sum over Matsubara frequencies in Eq. ( 12 ) and consider 
the zero temperature limit. We find 



FBZ ^ 



U 



1, 



4c/(p-|rg 



(13) 



which, together with the equations for the Qs, completes 
the set of self-consistency equations for the Qs and p for 
the rotor-disordered (MI) phase. For the rotor-condensed 
(SM/BI) phase, as in the case of a normal Bose- Einstein 
condensate, a macroscopic fraction, namely hq per lattice 
site, of the particles occupies the A; = state, which cor- 
responds to the lowest energy — iFglmax- The chemical 
potential (— p) is equal to this lowest energy. The corre- 
sponding equation becomes 



Therefore, there exists a critical value 
t'Q'f^2tQf 



(8) 



that separates a SL with gapless spinon excitations, 
around two distinct Dirac points at the Fermi level (for 
t'Q'j: < 2tQf), from a gapped SL phase (for t'Q'j: > 2tQf). 
The spinon dispersion calculated from the Hamiltonian 
T-lf is depicted in Fig. [3] for two characteristic sets of pa- 
rameters. 

To obtain the Mott transition line, one needs to self- 
consistcntly solve for Qg, Qg, Qf, and Q'y. This is achieved 
with the help of the imaginary-time action 



S = 



3 



iLjdrOj 



fja^rfja 



n 



(9) 



where /3 = l/^ksT) is the inverse temperature and kg 
the Boltzmann constant. We next rewrite ^ in terms of 
the fields Xj = exp(i0j) and constrain the normalization 
|X,(t)P = 1 via a Lagrange multiplier p. The resulting 
Green's functions for the fields X and / are 



G~\k,iv^) = -[j,l/u + p-\T^g\), (10a) 

G-f\k,lU3,,) = l^n-\T% (10b) 

where w„ and Vn are the fermionic and bosonic Matsub- 
ara frequencies, respectively. A change of ?7 — > J7/2 is 



performed in equation (10a), in order to preserve the cor- 
rect atomic limit [21]. Based on the form of the rotor 
Green's function, one may define the charge gap as 



(11) 



2\U{p-\Tl 



where iFaLax is the maximum of iFfll over the FBZ. The 



normalization of the X-field yields the equation 21 
(Pk 

'^Geik.ii'n) 



FBZ 



-1. 



(12) 



We consider the rotor-disordered and rotor-condensed 
phases separately. In the former case, we perform the 



2no 



(Pk 



U 



FBZ 



1, 



ir^l) 



(14) 



which, instead of Eq. (13), completes the set of self- 
consistent equations for the Qs and uq for the rotor- 
condensed (SM/BI) phase. 

In the rotor-disordered (MI) phase, one has Ag > 
{p > iFglinax) and no = 0, while in the rotor-condensed 
(SM/BI) phase, one has no > and Ag = 0. A second- 
order transition line is defined by Ag = and no = 0. As 
we show below, part of the phase transition curve that we 
obtained is of first order, indicated by a jump of Ag or no 
from some finite values to zero. For the rotor-disordered 
(MI) phase, one finds three distinct solutions with different 
types of spinon excitations: 

The chain solution has vanishing mean-field parameters 
on the horizontal bonds, i.e. = Q'g = 0, and thus 
describes a system composed of decoupled vertical zig-zag 
chains [corresponding to the bonds with a hopping t in Fig. 
[2]^a)]. A self-consistent calculation yields Qo = — 2/7r « 
—0.637, while Qf and p are functions of U. The associated 
SL phase in this case is gapless, since t'Q'j = < 2tQf. 

The dimer solution has Q f = Qe = 0, and hence de- 
scribes a system composed of decoupled dimers on the 
horizontal bonds of the honeycomb lattice [see Fig. [2]Ja)]. 
Self-consistency requires that Q'^ — 0.5 and Q'g = — 1. 
One may easily verify that Ag = U for this solution, inde- 
pendent of t' . This relation also fixes the value of p. This 
SL is gapped in the spinon channel, since t'Q'j = 0.5t' > 
2tQf = 0. 

A third (more general) solution in the rotor-disordered 
SL phase is continuously connected to the rotor-condensed 
SM or BI phase, and it therefore separated from the latter 
by a second-order phase transition. This solution may be 
obtained by starting in the rotor-condensed phase with a 
sufficiently small value of U, where the mean-field param- 
eters are obtained numerically. These parameters serve 
as the starting point for the next calculation step (i.e. for 
U + e, with an incremental energy step e > 0). One thus 
obtains a continuous line of solution in parameter space. 
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Fig. 4: (Color online) Zoom on the phase diagram in the vicin- 
ity of t' = t. For the notation and meaning of the symbols, 
refer to Fig. [l] The dashed vertical line denotes the isotropic 
hmit t' = t. 



The condensate fraction uq decreases with increasing on- 
site repulsion U, and the rotor-disordered MI phase is ob- 
tained when no < 0. Furthermore, the charge gap Ag 
increases from zero continuously when U is increased, as 
one expects for a second-order phase transition. 

In order to identify which of the above-mentioned so- 
lutions is chosen by the system for a particular set of 
parameters {t' /t,U/t), we have performed a free-energy 
analysis of the three solutions. The free energy per lattice 
site reads £ = —{kBT/2Nc) InZ in terms of the partition 
function Z and the number of lattice sites 2iVc, and one 
has, in the zero-T limit, in which the free energy becomes 
the internal energy 



£^£f+£g- [2tQfQe + t'Q'jQ'g + p], 



where 



FBZ 



<Pk 

it' 



-/I 



is the free energy of the fermionic part, and 



FBZ 



2n 



(15a) 



(15b) 



(15c) 



that of the bosonic part. Note that in the condensed 
phase, p = iFglmax- The last term in equation (15a) ac- 



counts for the correction to the free energy due to the 
dynamically unimportant constants, which have been ig- 
nored so far. 

For a particular pair of {t'/t, U/t), the true ground state 
is the one with the smallest £. The transition curves re- 
sulting from this free energy analysis are plotted in Fig. [l] 
(lower thicker green and dotted line) , and Fig. |4] for the 
vicinity of the isotropic limit {t' = t). One finds that, for 
t' < 1.005t, the system undergoes a second-order phase 
transition (within the third solution) from the SM phase 
to a gapless SL phase (GLSL in Fig. [l]) upon increase of 
U/t; for 1.005t ^ t' < 1.019t, the system first experiences 
the above-mentioned second-order phase transition from 
a SM to the GLSL (third solution), and then a first-order 



phase transition to another gapless one of the chain solu- 
tion [GLSL (Chain)]; for 1.019i < f < IMt, the second- 
order phase transition disappears and the system jumps 
directly from a SM phase to the GLSL (Chain) via a first 
order phase transition. Finally, for t' > 1.31t, the system 
undergoes a first-order phase transition from the SM/BI 
phase to a gapped spin liquid of the dimer solution [GDSL 
(Dimer)]. The AF transition, which cuts the GLSL-GLSL 
(Chain) transition at t' ~ 1.005i, but which is otherwise 
well above the SM/BI-SL transitions described here, is not 
treated within the current slave-rotor description. In prin- 
ciple, this AF transition may be described within a slave- 
rotor theory by reintroducing spin-charge correlations into 



the Hamiltonians (5b) [23]. 

For the isotropic case t' — t, our result of the second- 
order transition point Uc/t ~ 1.68 agrees with Refs. 18 
[20) . In this limit, the free-energy analysis indicates that 
the gapless SL of the third solution is the lowest while the 
gapped SL of the dimer solution is the highest in energy 
among the three solutions. This is in contrast to the QMC 
calculations |16 , where a gapped SL was identified as the 
true ground state. However, only an upper bound of the 
free energy of the dimer solution is provided by our free en- 
ergy analysis. This is especially the case in the isotropic 
limit, where a discrimination among the three hopping 
parameters is unphysical. Indeed, possible kinetic dimer 
terms (e.g. resonating dimer moves around a hexagon) 
could lower the energy of the gapped SL of the dimer 
At a general value of t', free energies are 
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modified by these terms such that they are adiabatically 
connected to the (also modified) isotropic case. Similarly, 
the gapless chain solution may be unstable to an infinites- 
imal interchain coupling that, in the context of the square 
lattice, is known to open a spin gap f25'. An analysis of 
kinetic dimer terms and the stability of the chain solution 
is, however, beyond the scope of this paper. 

The predominating SL phase in the limit t' ^ t is de- 
scribed by the dimer solution, where both Qf and Qx are 
zero. This in effect is equivalent to setting t — 0, such that 
only sites connected by horizontal bonds are coupled by 
tunneling (Fig. [2]) and the whole system reduces to a set 
of decoupled dimers. It is therefore sufficient to solve the 
Hubbard Hamiltonian for two interacting fermions occu- 
pying two lattice sites coupled by the hopping amplitude 
t' . The two-particle Hilbert space obviously reduces to six 
states as the Pauli principle forbids fermions with identi- 
cal spin to occupy a same site. By the same token, states 
where fermions with identical spin occupy different sites 
are obviously zero-energy eigenstates as hopping is Pauli- 
blocked. The relevant two-particle Hilbert space is thus 
spanned by | t,l), I tl, 0), |0,tl), and | |,t), where the 
first entry describes the occupancy of the A site and the 
second one that of the B site. The diagonalization of the 
Hubbard Hamitonian then yields the ground state 



sm ( 



\S) 



cos ( 



n,o) + |o,n) 

V2 



(16) 
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where \S) = (| t,i) - I l,t))/V2 is the dimer smglet and 
(f) — (2a + 7r)/4 with tana = U /At' . The first excited state 
is the triplet state |T) = (| i) + I l,t))/V2 and is sepa- 
rated from IV'gs) by the gap Agx = [^/iP + W? - U]/2. 
The ground state is therefore essentiahy the singlet state 
|iS), with an admixture of states with double occupancy 
the weight of which {coscf}) vanishes in the laige-U/t' 
limit, where the gap is dominated by the exchange energy 
AgT ~ At'"^ /U . Because IV'gs) remains the ground state 
over the whole range of values of U /t' , there is just one sin- 
gle phase that continuously connects the BI to the GDSL 
phase. Notice that, within the mean-field U(l) slave-rotor 
treatment, there is not such a continuous connection be- 
cause a can only take the values (for the BI) and 7r/2 (for 
the GDSL), double occupancy of a single site being ruled 
out in the MI phase. Furthermore, only at the asymptotic 
point, U = oo, do the singlet and the triplet states become 
degenerate, and the AF state, which requires a degeneracy 
(and a superposition) of \S) and |T), may be formed. In 
the intermediate region, the gap Agx protects the ground 
state IV'gs) which is only marginally perturbed by a small 
inter-dimer coupling mediated by i ^ t' . Within this 
simple dimer picture, one may therefore understand the 
gapped SL [GDSL (Dimer)] in Fig. [I] as being adiabati- 
cally connected to the state IV'gs)- 

In conclusion, we have investigated the repulsive 
fermionic Hubbard model on the anisotropic honeycomb 
lattice. This system could be experimentally realized by 
loading ultracold fermions at half-filling in a honeycomb 
optical lattice. Beside the SM, BI, and AF phases, which 
are readily obtained at the mean-field level, we have used 
a slave-rotor description to show the emergence of var- 
ious SL phases. Two gapless SL phases may be found 
from a self-consistent mean-field treatment of the slave- 
rotor theory: a chain solution that consists of essentially 
decoupled zig-zag chains and a phase that is connected 
continuously to the small-C/ SM/BI phase. The latter is 
found in an intermediate coupling regime 1.68 < U /t < 2.2 
for the isotropic case, t' = t. Most saliently, a gapped 
dimer SL phase may be stabilized at rather low hopping 
anisotropics, t' > 1.31, and is adiabatically connected to 
the large-t'/i limit, where the system may be described 
in terms of essentially decoupled dimer states. One may 
speculate that the dimer picture, if modified by an ad- 
ditional kinetic term that allows for dimer moves, yields 
insight also into the physical properties of the system in 
the isotropic limit t' = t. 



We thank A. H. Castro Neto, F. Crepin, B.-G. Englert, 
J.-N. Fuchs, M. Gabay, N. Laflorencie, K. Le Hur, G. 
Montambaux, F. Piechon, and M. Rozenberg for fruitful 
discussions. We acknowledge support from the France- 
Singapore Merlion program (CNOUS grant 200960 and 
FermiCold 2.01.09) and the CNRS PICS 4159 (France). 
The Centre for Quantum Technologies is a Research Cen- 



tre of Excellence funded by the Ministry of Education and 
the National Research Foundation of Singapore. 



REFERENCES 



A. H. Castro Neto, F. Guinea, N. M. R. Peres, 
K. S. Novoselov, and A. K. Geim, Rev. Mod. Phys. 81, 
109 (2009). 

M. O. Goerbig, (preprint) arXiv: 1004.3396 

V. N. Kotov, B. Uchoa, V. M. Pereira, A. H. Castro Neto, 



[2] 
[3] 

[4] 

[5] 

[6] 

[7] 

[8] 

[9] 

[10] 

[11] 
[12] 

[13] 

[14] 
[15] 

[16] 



[18] 
[19] 
[20] 
[211 



[22] 

[23] 
f24" 



F. Guinea, (preprint) |arXiv: 1012.3484 

J. E. Drut and T. A. Lahde, Phys. Rev. Lett. 102, 026802 

(2009); Phys. Rev. B 79, 165425 (2009). 

X. Du, I. Skachko, A. Barker, Eva Y. Andrei, Nature Nan- 

otechnology 3, 491-495 (2008). 

J. Gonzalez, F. Guinea, and M. A. H. Vozmediano, 
Phys. Rev. B 59, R2474 (1999). 

I. F. Herbut, V. Juricic, and B. Roy, Phys. Rev. B 79, 
085116 (2009). 

S.-L. Zhu, B. Wang, and L.-M. Duan, Phys. Rev. Lett. 98, 
260402 (2007). 

K. L. Lee, B. Gremaud, R. Han, B.-G. Englert, and 
Ch. Miniatura, Phys. Rev. A 80, 043411 (2009). 

T. Bourdel, L. Khaykovich, J. CubizoUes, J. Zhang, 
F. Chevy, M. Teichmann, L. Tarruell, S. J. J. M. F. Kokkel- 
mans, and C. Salomon, Phys. Rev. Lett. 93, 050401 (2004). 

P. Soltan-Panahi, J. Struck, P. Hauke, A. Bick, 
W. Plenkers, G. Meineke, C. Becker, P. Wind- 
passinger, M. Lewenstein, and K. Sengstock, Nat. Phys., 
doi:10.1038/nphysl916. 

C. Lee, X. Wei, J. K. Kysar, and J. Hone, Science 321, 385 
(2008); M. O. Goerbig, J.-N. Fuchs, G. Montambaux, and 
F. Piechon, Phys. Rev. B 78, 045415 (2008); V. M. Pereira, 
A. H. Castro Neto, and N. M. R. Peres, Phys. Rev. B 80, 
045401 (2009). 

P. Dietl, F. Piechon, and G. Montambaux, Phys. Rev. Lett 
98, 236405 (2008); B. Wunsch, F. Guinea, and F. Sols, New 
J. Phys. 10 103027 (2008). 

S. Sorella and E. Tosatti, Europhys. Lett. 19, 699 (1992). 

N. M. R. Peres, M. A. N. Araiijo, and Daniel Bozi, 
Phys. Rev. B 70, 195122 (2004). 

Z. Y. Meng, T. C. Lang, S. Wessel, F. F. Assaad, and 
A. Muramatsu, Nature 464, 847 (2010). 

S. A. Jafari, Eur. Phys. J. B 68, 537 (2009); W. Wu, 
Y.-H. Chen, H.-S. Tao, N.-H. Tong and W.-M. Liu, 
arXiv:1005. 2043^^2. 

S.-S. Lee and P. A. Lee, Phys. Rev. Lett 95, 036403 (2005). 

M. Hermele, Phys. Rev. B 76, 035125 (2007). 

S. Rachel and K. Le Hur, Phys. Rev. B 82, 075106 (2010). 

S. Florens and A. Georges, Phys. Rev. B 66, 165111 
(2002), ibid. 70, 035114 (2004); for a pedagogical re- 
view of this technique, see E. Zhao and A. Paramekanti, 
Phys. Rev. B 76, 195101 (2007). 

G. Montambaux, F. Piechon, J.-N. Fuchs, and M. O. Go- 
erbig, Phys. Rev. B 80, 153412 (2009); Europhys. J. B, 72, 
509 (2009). 

See e.g. E. Zhao and A. Paramekanti, Phys. Rev. B 76, 
195101 (2007). 

R. Moessner, S. Sondhi, and P. Chandra, Phys. Rev. B 
64, 144416 (2001). 



p-6 



Emergent Spin Liquids in the Hubbard Model on the Anisotropic Honeycomb Lattice 



[25] For a review, see I. Affleck, J. Phys.:Condens. Matt. 1, 
3047 (1989). 



p-7 



